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SCm^TlC R8H?ES©IThTIC® or DSSXOS process 



SELECT PRINCIPAL PARAMETERS 



DETERMINE COORD GcATES OP 
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<• a' jt h £t< bo an .J. to % j< , •... _»• c ,a >. .... 

rffer ro ■— s design©.*. gract** s’jou^h o 'h uty ^ , , u . -y =. • l y. a t 
*w.ipaL4.: . .Pith efficient coo ma/; usage . j<te .•iif.vr . ly «*ir,t tist 
lu^ern*! pperatica of the pro^stui xe aou sectssary . cut t, acwuilai 1 . 
he ca*a: the user will gali caxi loanee .to the imre-ei. operation. 

Eli® seem* of plating ore Dum-esred in age t.adiag ordc ' :-.:co: 

Jha heel plate center up the shell aids; and contiiauxig frcm. the 1st 
deck to the inner bottom, if thr.t exists; and finally laroucji "he 
bulkhead seams. Tbe shir has keen broken down in so see 'loi.-,-* end the 



longitudinal* within those sections ar c - indexed in cue fo.-lowiag 
'■aanaer. 



The bottom she. I is xahelc-o -3 section f 1 - uic first 
longitudinal being indexed as {10.1} u^d contix.al.-a tier the total 
number of losgituditais for taat sectlja. The aes: 1 - s«<»tj,.n up tne 
shell side is labeled section #11, the next is #12. ' or»tiii'.siug until 
the a.5 in deck is ro&eked. The desks are jab*': Jed as aiolion #1 lor 
the iredn deck, section #2 for the second, continuing fox 11 the deois, 
the inner bottom is indexed as axe more than the nuatwr of decas .> 



Again the longitudinals fellow the 3 esr: xogie as two fils, Therefore, 
(3,7) signifies the seventh longitudinal on the tbi.-d deea 'The 
bulkheads begin with the Index of ‘2J“ fn the "Irut bulkhead, and 
continue for all bulkheads. Again the Icrgitudir.als ere labeled ae 
before. Before operating the program, a sketch sheuid he prepared 
so any confusion will bo a? deviated. 
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The actual operatic^. begins wito. the loading of the 
urogram into the computer core aa shewn belou. 

LOADC-C !<H STF PIS TA? £3JB CYG BUK SSi'B 
S’j^t 2- 

Presently the jrograifl is set cc &ei>*g& c :.erefca.at snip, and 
*11 principal dimensions are arranged to fit A.E.S., chips , Tks Input 
dace existing now cen be modified. Thseaiaguter vili print a causHSit 
shoving what symbols .can be changed at each interval. and the changes 
Sky be sna&e. 

This is dcno by setting toe symbol representing tie quantity 
to 'be clanged equal to t ae .new value. For eraEtple, (UP * 500) aieans 
that the length between perpendiculars will be 5 GO rt„ 

The information telling which sections are to "oe transversely 
fraued is contained in a vector called TrKAKSV » If sections 1> 2, and 
11 are to be transversely framed, the array is filled by typing 
l! TRARSV(l) - 1, 2, 11”- Any section nuasber may be included on the right 
side of the equal sign. 

Longitudinal bulkhead information is contained In & two- 
dimensional array called ,: SHD (l, j)”. The rows (I) are the sequence 
of bulkhead (i.e. , firsts second) while the columns (j), contain 
data describing the location of the bulkhead. If the operator wanted 
two bulkheads , and the second were 10 ft fran the center lias , lying 
between decks #1 and $ 3 , he would type "3133(2,1) - 10, 1, 3.'* 



k 



An inner bottom is provides by typing; IB - $YES$; if no 
inner bottom is desired then type t IB - $H0$. As each interval of 
data inpit is finished, the operator must type an asterisk (*). 

Befei'encs to the sample session at the consols will clear up any 
confusion found in the scanner in which date. Is introduced. 

at two internals la the design of the coordinate system, the 
operator may atop ths program to look at the scantlings. This provision 
is only to be used if the coordinate system is being revised, and 
provides a quick look at the locations . In any event, they wil l be 
printed out at the end of the program, so it is best to skip them now. 
Control is removed from the operator until the end of each 
cycle, whereupon he can .make his decision as to which set of thicknesses 
to use for the next cycle; or if the design ia satisfactory, fee may 
choose the final thicknesses and stop. The final scantling design is 
printed end the program is ended. 
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SAMPL? 



,!O l ! AJ. OlHSOlP 



londno nn pit stf tap eye buk oru! sub 
2021 „ 7 
rc!ii i r*n . 



THE VALUES BELOW MAY BP ALTERED OOi 



LRP, R, P, H 1 /! r >KS,HfEL,n p 'R| SF, FOP P, WEBS P, S ! '‘'MAY , F , fiAWr'f; 
YRAS I 0, DBAS » 0, AR/vS ! C, SBA3 I 0 , OX,. OB, ! P ,TRA'I3''( ] ) 



lbp*550,h»75,d*39.3,ha27.5,cx=.98, eb=.67,f rsp=30,wnbs ->=3 0.iv«ks- 
nbhds=0, i b=$yos $, 
transvCl) =1,2 ,3,4,10,11,1 2* 



THE VALUES BELOW HAY RE ALT EREP HOW 



S I 0 LI M , RM , R, H ATCHO , C VK , 7 R F. o 



cvk=3.93* 

EOR ST I FFEHFRS,TYPE l-HOMR,T v PF 0 



0 

FOR SEA !S, TYPE 1, HOME, TYPE 0 
0 
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COORD 3 Y COORD SCANTLINGS SECTION MODULI* 
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DECK WO : 5 ! S TRANSVERSELY PRAMF 



OF SCR? P T 10 H OF IWTFRMA L ') P"RAT l • 
FLO j CHARTS 



*» 16 *» 



F'l W 3' 



■ v . 




A decision branch 
S’ "vheneves^' sta cement) 




NO 




An iteration; for ex£mpj.e, the ’through loop’ 
’■•/here I is initially set to 1 and tin itera- 
tion is performed until ] is prsater than 
TOTS’fP. » 



AREA 

CYCJS 

METER 



Tfee enclosed raluea are calculated. 




Subroutine (BUK) is executed. 
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<tr> 



Continuation from previous page 



Coat isuat. ion co next page 



Print out of data 



Read tn Eata 



Beginning or end of each subroutine 



=x* 
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(l?8. ) a t a*c cu.ii.tog program tt • executes <fae 

•nubroutines and soHU&iD.s all of the necessary input data for the 
< pera r ion of che program. 

(MN . ) first clear': all ,3 insulated tapes o i any existing data 
comput ad from previous runs and then set? the values for the folrowingr 

Young's reodilri, length let, veer perpendiculars, beau, iepth, 
draft, number of decks, heel angle, deu&rise, frtuae spicing, web 
spacing, midship section coefficient, height of center vertical keel, 
block coefficient, yield stress, number of by lift cads, the depth, area, 
height of neutral axis, moment of inertia and section modulus of the 
ba3lc stiffener, bending moment , limiting bending stress, bilge radius, 
number of longitudinal bulkheads :ai the midship half section, and the 
width of cargo hatch. 

Tire required ship section modulus io determined according to 
the Load Line Regulations and is calculated by the subroutine (PLQSEC. ) 

To change ary of the above existing data, there are two program 
stops which permit new .input data to read into the computer. Reier 
to the section on program operation for further clarification cf the 
main program. 
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J 



M 



ABASIC 

3 

BILL 

BM 

CB 

CBILSE 

CBR 

CVK 

CX 

D 

DBASIC 

DDRISE 

E 

EM 



FLBOT 

FRSP 

H 

HATCHO 

HB 

HSEL 

IB 

XBA3IC 



c~ j 0 'ALL AHID I'lll 

Array of numbers from (I), Sect-. 43,15 ~ 17 (c j (5 ) 
used in calculating ZB3Q. 

Area of basic stiffener. 

Breadth of ship. 

Linear length of bilge. 

Bending snoment of ship on wave , 

Block coefficient. 

Height of bilge center above keel. 

Standard carrier 
Height of inner better. 

Midship coefficient. 

Depth of ship. 

Depth of basic stiffener. 

Dead rise of ship hull. 

Young ‘3 modulus. 

Function used in a buckling formula later in the 
program. 

Length of flat portion of bottom shell. 

Nominal longitudinal spacing . 

Draft of ship. 

Half breadth of hatch opening. 

Half breadth of ship. 

Heel of ship used in calculating incline'! •water head 
Inner bottom 

Inertia of basic stiffener 
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LIST 0? SYMBOLS C fiLCfOI ATED Tk b'^i 
“* £ 


LBP 


Length between pexpendiciileirs . 


.LBEDS 


Ihasber o;:' longitudinal bulkheads on aall section. 


HDX3 


NuEiber of sti*orgth decks ol ship. 


PHI 


Angle of deadriae. 


j 


3*1416 


fi 


Bilge radius. 


sigum 


Limiting primary bending sts-esa. 


SI©IAY 


Yield strength of steel used. 


TRANSV(l) 


Array consisting of the numbers of sections to be 
transversely framed. 


TWD 


'Tween deck height. 


WEBSF 


Spacing of web frames. 


YBA3IC 


Distance fro® shell to aeutial axis of bs.sic 
stiffener. 


ZREQ 


Section modulus of ship required by (1). 
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Til?. 



(TAP.) is a subroutine wU£c . reads and writes PEOPCO 
and XLL(l, J) onto simulated tapes i*t various points in the program in 
order to conserve storage locations In the memory osa of the computer. 



(REQSEC. ) is an internal function which calculates the 
required midship section Modulus according to the Load Line Regulations . 
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£EAL\ 



{EKA£o ) is an interns.! function which. e&leul&tea water head, 
using the horizontal (ir.il vertical G-sordjinstss of s. point oe the bull 

geonetry as inputs. 

Thf internal function is executed by (STF. ) t&d {.PLT. }, and 
calculates two values cf bend for each coordinate; fit* 1 head when the 
ship is heeled thirty degrees and the head in way of the wave erset. 
The larger of the two values is selected as the water head and it 5 a 
not to he less than 0>>£5 feet. 

LIST OF SYMBOL? CALCUIATES XS (ffitfj ) 

33iCI£J> Water head on a structural iroscdJer when the at .ip is heeled 
30 degrees. 

WHD Water head on a structural mecher in way of s>. wave urest 

HM1S Mininsua water head cf ?.2|5 fest. 
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Sis! t e rkikl 



f'UKiCT \«}M 



heap,(x,y') 



e 7~X A RT 
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«p 

„L 



INCLHfe -(tt-Y* 

V/HD *$>* LEP.EO.& -V 
HMiN =. 5v 1.7- 
£X~- HMIM 



;;«Ces„(Ht£L) 



V 






= W Hx^> 




, X 

1M CL 

»Lj 1 #*- — — r 



NO 



HiA/\x - 
5NCLHO 



j 



v 



EH1> 



- S7 - 



SXP. 



f iiTJ o } csil cult* t»es the IfiSJgitudlsiiil frssru k pac „*& each 
sect im ef the ship.* determine® the horLtontal and vertical cfmr&tnateg, 
iss&awsred fxvsa the keel, of the laagitudiwal stiffeners, aad calculates 
th® water head an each stiffened loootisa* 

The girths of e®.ch section are first computed aa& then a 
nesaiial longitudinal frame spacing 3 s detemined so that the longitudinal 
©tlffesxerg tore eqv.slly spaced an each section sad as near as possible, 
but not exceeding, the input frame aptclrg. 

Hext, the longitudinals ©re located cn e&oh sectiea based cn 
the spacing **v that section, ©ad the XII. {horizontal ; and m. (vertical) 
coordinates are calculated. With these coordinates, the water heads on 
each stiffener are computed by the subroutine (HS&'Q, )* 

Fvr sections that are tmisv&roely fraaisd, thus having uo 
longitudinals, XLL, YLL end the water heads for tile sections are oet 
ecpial to zero at the end of this ©ubroutine. 
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..YbBOLS CAICOIAXKD 


ALPHA 


Angle measuring distance of inner bottom fro a center 
point of bilge turn. 


B'JIST 


Incremented son; of frame spaces oo ewe bulkheads. 


BILGrCH 


Incremented sum of fiams spaces on bul^ft turn. 


CC)UKT 


The accumulative sum of longitudinals cu shell side. 


D GUTA 


Incremented angle between longitudinal upper bilge 

turn. 




Array of the deck heights above the keel. 


FEDlV 


The first estnsate of the number of lougitudinels 

in a section. 


GAMMA(l-J) 


Angle of longitudinals on bilge turn up tc inr.er 

bottom. 


GIRSH(I) 


Array of all the section lengths. 


HMAX 


Maximum water head chosen. 


HM2H 


Minimum water head possible. 


IBHB 


Inner bottom half breadth 


IHCLHD 


Head on longitudinals when ship is inclined. 


I£<pjT(l) 


Astray counting the longitudinals in each section. 


LHD(I,J) 


The water head on each longitudinal. 


LL10 


Accumulated sum of the longitudinals in Section 10. 


im 


Accumulated sum of the longitudinals in Section 11. 


LLDK 


Accumulated su a of the longitudinals on the dec Is . 


LLLB 


Accumulated sun of the lergitudinale on the inner bottom 
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SYMBOLS JAICUMTSS IN/STd.) 
2 «» 


S(I) 


Array of all the frame spaeings in each section. 


SJ5IST 


The incremented frame spacing of each section on 
the shell side. 


SDVAK 


An integer which has the value 10X3 when the ship has 
no inner bottom and the value 10X3+2. when an inner 
■bottom exists. 


SIDE 


An incremented length at bilge hotted* similar to 
BIJjSTHc 


swam 


An incremented length in section 11 on the bilge turn 


SUM 


A quaait.v describing several incremented totals used 
through out the program. 


TDIST 


Incremented frame spacing seen on the top deck. 


TEMPL 


The distance between the last longitudinal on the 
bottom shell and FLBOT, 


t^tall 


The total number of longitudinals on the half section 


WHD 


The wave head on the longitudinals „ 


XLL(lgJ) 


"X" coordinate of longitudinals . 


YLL(I,J) 


"Y" coordinate of longitudinals. 
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fpxic. ) 



is a ff;«brc it^.as which deiormiuea a,.l gc-csnetiic 
properties that depend upon the plate aeons. When execution of this 
subroutine has been ceEjleted. 'idle ci> v rdinate system iot all. longi- 
tudinal fcatiauous material 1- -«he *n ids hip section is ; Lxed atA the 
program Jjs ready to beg.m tie scant! tag selection procedure . 

The first step ii ihe execution of the t! PLT is the fixing 
of plate widths. The keel, center, margin, sheer, toj and stringer 
stroke widths have been calculated in Ref, ( j according to iS5 rule.:. 
The relations obtained art! used by tills subroutine to x>rese\; the values 
of these plates. Hexfc, a standard plate width based upon ship length 
is selected for all remo.is.iixg plating. 

When a standard else has been selected, the next step is to 
test all seam location* for interference with the stiffener; vnose 
coordinates have been fixed, Ihe keel, ^late is checked first and 
corrected if interference occurs, then the center Strabo is made to 
equal the keel width. The remaining standard plate widths must now fee 
changed to a nominal size which yields an integer* number of plates on 
the shell, the decks ana on the bulkheads. Thus, every otrc.ke width 
has been determined. It remains now to check each one for inter** 
ference with stiffeners end tc store away all data a&edea by subse- 
quent subroutines , 

The shell, seams are first examined. The distance dong the 
shell from the center of the keel for each s brake top is found: checked 
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against ;=tiffeae.i.- in-arxcrc.K:^. end ca^jed. :'.f scc>,&s^jfy. Eext tac 
coordinates of the s earns an:>. titj vx ter he«.£ o.- then. is calculated. 

The area, moment and inertia constants irs also calculated for use 
later in the program « 

Deck, inner bottom, and bulkhead plates axe t.Miaed In a 
sequence similar to that of tie shell. -The disianc of each seam from 
the deck side or bhd. bottom is found and checked agalr.st interference. 
Then the coordinates, water head, and area moment aai inertia constants 
arc found. Once these quantities ai*e calculated, the program returns 
control to the main program. Before this is done, however, it is 
possible to obtain a print out of the data calculated in this 
subroutine. Refer to the chapter on program execution for the correct 
method. 
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J .0 0 -- . V** 1 - 1'Cj.i 'z - . ■ U J . (; C - .11 . L'j! } 


AC ,J) 


Area -om ’teat or platex .ici.nnVet. 

length of ? lets. 


BDSTP.(l) 


Number of t teak as on bill. 


COIiTl 


Key to prix.t additional clc ia 


CTRSTR 


Width of center sc rake on inner bettc* . 


FKSTR(l) 


Number of ete-akes on I"th deck. 


Fift(j) 


ingle between seam- on bilge, sirs lor to " SAidt in ’s^F 1, , 


HSEAM(J) 


Water head on ' J i! fch ssem 


XC(J) 


Inert-La cec3taat. of plate 


keel 


feel plate width 


Margin 


Width of margin plate on inner hot tom 


MC(J) 


Moment constant of 'J’th plate 


REMP.UT 


length of all plating on shell except iCELu, 3HL2R, T? SIL. 


SEAM(J) 


ifee linear distance from the top seam of tee V ' plate 

to the Keel . 


SHEER 


Width of sheer plate. 


SHELL 


Linear length of ail plating on shell 


SHTjSTP. 


Number of plates on shell. 


3PANST 


The plate to be associated .d.th each stiff* aer. 


STDPLT 


A nrninal assumed standard width of plats to fee- dependent 
on 3hip length* 


stkstr 


Stringer plate 


SUMl 


Accumulated sura of GU333l(l} v.p tae side freer the keel. 


TODK 


Number of plates on ail the decks 


TOTSTR 


Accumulated sum of ail the plates on the half section. 
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TPSTR 



WBDSTRCl ■ 
W0KS1‘K( I ) 
WSSIST 
X SEAM (J) 
Y SEAM (J) 



oi l-s'Ilco d 



■ . -'V - i : 

T»«.. iSMW t 



‘Yop atrake 

ISaainal width, o: each plate ‘ I : tb oh*?.- 
Nominal width of* each plate cm the • I -to c.cc-1. 
Hoadnai width j? each plate cn 3 hell. 

"X” coordinate af ' cenn 
"Y" coordinate ef bo,u.i 



✓ 
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NP 



AlOSBF. 



{AM1STP. ) calculat s# tta total euv'S.., uasimt end inertia 
of the stiffeners at oat the base .line. 

Using the value oi prep ore i coal 1 fcy constant for each 
stiffener location detwiaiaccL ia {STPHER. ), the arc®., maaent and 
inertia of each stiffener ie calculated. The result las value* fox 
each stiffener can be added, ssSnce they were tahea about the s&ae 
bane line, thus giving; the total sure®, aesaent end inertia of all. the 
stiffeners , 



LIST OF SgjjgOLS CALC mA'Effi IH (AKlSgg, ) 



ASEF Total cross-sectional area of the lougifudiml stiffeners 

at tha midship section, 

2STF Maaent of inertia of the cross -aectioafcl area of the 

laagltudiaal stiffeners about the base line. 

MSTF Internal bending accent of the ercoe-secticaal area of the 

longitudinal stiffeners about the base line. 



o 
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AMIS' 



5 , 
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AMXPI2' 



(AM2?L£ „ ) calovi.«,te a the total area, secaient and inertia of 
the plates about the base lies. 

Using the values of area cons taut, consent constant and inertia 
constant calculated in (P13?.), and the current value of plating 
thicknesaes in the eub= cycle, the area, consent and i»e?rtia. of es’.oh plate 
about the be.ee line is calculated. The resulting values arc added, 
giving the total area, docent tat> inertia of r_U plates about the 
base line . 



USX OF SYMBOLS CAICUIATEB SH >AMI?LT. j 

APLT Total cross -nee tion&I area of the plat leg at the isidship 

section. 

IPIT Mecjent of Inertia of the eross^Gectloaal area of the 

plating about the bane line. 

MPLT Internal or resisting bending ccssent of the or®ss"Sectiona.l 

area of the plating ebout the base line. 



AM 1 PlT. 



(M’TART j 



J_ 

5 A per 
Mper 
XPCT 



X*i 



j 



APLT 

mpltI 

TFirrJ 



EWD 



3 
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JJ • — .! . : 

J-Oi Xt’ji . , _ «. • . 

tiie fc^birutHae v . •* ,' , -tL -j. 

L.J Ldj~ Lz.yz.: v. ' c . 4 ... . ... 1 

proportiomaity eat»ts*ri c: tfee b$u»i5 :•• y _. 
required section or « platc-rf 5 '- r.: er sou- . •£.< : 

LIST 0.1‘ il A T;Q SIS (in'jffTT- j 

TS4? IXisny variable equal to &he vc^pc*rtionc ..Iq 



. 1 . 

3 .. . > la 

^o.. . ; ..> *. 
cr. . 



• constant. 
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w» .C • * CA .L'.VtSb 

'izt f Jgitctlir ..i *..■:. ie r '.q’.:u, . c u • .i. c- .. 

f ion to '''iiket&jQi the i -xxte^m ntjsst. • 



^Leaovev (STvfciai.) iu extcuk-J: i>y {CTC. ) . the ecr^&iK&csa v 
the stiffeners, water £*&ds, plate tkie.Vrieijs-v-i arseo^i-' o. st* vita each 
stiffener location, and the stresses f<:a the current sub-eye !•;• arc 'rat- .. u 
It 5,s the purpose of this sabroutii c to determine toe repaired stiffen * 

scantlings « 



A scale factor os- proper locality constant, liliPOO ; based oa 
the geometry of the basic stiffener end an effective plate ailth of 
60 x thickness or the longitudinal frawj spacing, whichever is less, la 
used to determine the geometric properties of the required stiffener. 

For each stiffener .tocatica and its associate -i plate 
thickness, (STFHOER.) calculates fifteen values of section modulus for 
each plate- stiffener combination, by using fifteen different values of 
PROPCO. In effect, a curve of section modulus versus proportionality 
constant is generated for esch stiffener location. The actual section 
modulus required to resist the stresses in the- current sub-cycle is a 
point on the curve. To design a ship that approximately meets the 
American Bureau of Shipping Rules , a section modulus required uy the 
Rules, MODABS, is Initially selected as the ufxxi&m required section 
modulus of the platen stiffener combination. With this initial value- of 
MODABS, the subroutine (ISTSRP. ) interpolates the section modulus versus 
proportionality constant curve to determine the PRO.PCO associated with 
the initially selected value of MODABS. With, this value of FROFCO, 
(STFSER.) calculates the properties of the stiffener-plate combination 
for that stiffener location. 

«r> 
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J. — .. 4 . 5 i,. < v t ") ~j - f ■' ' / \. 

tfct *-» 'kiv.'U. x 3 ^** vx * k.ia« .nt. j f j » r c c * i v, j^c *$ 

~o‘ nula iui cc:io\nc»!l ^ • •*■■ _ .1 u..i-_:_ 1 . ...A ju c.ii . *j uou j. 

Sie section modu-Lua, mil, _i ou Tiruoshcr -j'q fosstula : -• n funct’.oc. c 

the geca&ctric properties 01 tb- tl:vU -stiffener coubicntiou and fu 
applied loads. Thus, corresponding to the iir.iW.al .•a.'ue of l O&'ci-i ard 
PROPCO , a valiie of section modulus. 22 IM, is oaiciuitt£ by Tlaoslienbc ■ -3 

Formula . Three possible situations exist at this point: 

I. ZT1M equals MODABS 
II o 220M less than MODABS 
III. ZT1M greater than MODABS 

In case I, when 2T1M and MODABS differ by not more than C.J5 > 
MODABS is considered to equal ZT1M, and is the required section modulus 
for the pia t e « s t if fenc-r conbination. 

In case II, since ZT1M is less than the ABS value, MODABS 
governs and the value of PROPCO associated with MODABS is the required 
one. 

In Case III, when ZH.M is greater than MODABS, ZT1M isust be 
the governing section modulus. But the value of PROPCO i®s determined 
from the initial value of MODABS and does not correspond to ZP1M. A 
new value of MODABS. is chosen such that the value of PROPCO determined 
by this new value of MODABS also corresponds to the value of ZP1M. To 
converge on a solution, the mean value between ZT1M and MODABS is 
selected as the new initial value of MODABS. A greater value of PROPCO 
results, and the two section moduli, 2T1M and MODABS based on this new 
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6l 



-4- .»-• iu' 



TLiu ubroutin^ ii. c.'LlzL d„ CjlV..> Ij. vt.'wi o rcitae th* 
thieknuuc 2.2 c... iatjiiq voquLrs . •••<, ,v .si .c ./ xer l.e>^c. 

2^.ch plate has two possible gcccr-tciC pm>_ ‘eaen •: i«,.e ; • . - o 

•the transversely frame* end loagitudiaally framed !^:iionj each s-e&.iw 
a separate stress distribution. She subroutine (M/-_'C3Tr.) ct-lled tw 
calculate the ssaxiir.in stress associated with each r„z\i on- -ulus wil_ 
be the stress upon the plates in the transversely iTccaesd Etctxons. 
Otherwise, the stresses calculated in (STBESl.) are ucad as a basis 
for the longitudinally framed sections. Those stresses ix (STBESl.) 
are associated with each seam. The stress on each plate is deter* 
mined from the maximum of the two seam stresses. 

The array CCk is next filled with the planing thicknesses 
based upon formulae that approximate ABS rule requires entr (Kef. Ho« )„ 
using the stresses described above. These thicknesses are those 
necessary to meet the stress schedule set up in the current sub*cycle 
being run in subroutine (CYC. ). 

Several different x r alues are used to calculate the 
thicknesses of each plate. Since the program is currently set to meet 
merchant, criteria, some plate scantlings depend upon the location, and 
upon the ABS rule requirements. The formulas for the center plate, 
keel plate, sheerstrake, stringer plate, and top strako are presently 
all special cases designed to meet AB-S requirements. A re-orienty.tioa 
of the long side in the longitudinally framed section requires a new 
fonnula for calculating CO^ This is triggered by the ore gram when a 
section having no longitudinals is encountered . 
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J 



T... - (CO. + x T )/ (WT - 1) 

A 



<- " c - - *■* 1" tdruc i »l*'u bated j» the last sub-cycle 

in subrc-.t ■ e ;r . i- 

1 1 ) 

T^g = Hew 'thickness of ith plate. 

T.^ = Old thickness of ith plate. 

CO^ - Thickness of ith plate just calculated based 
upon local lateral load. 

WT = Weighting function associated with ith plate. 

The weighting function is used to prevent the change in 
thickness during each sub-cycle frosi becoming too rapid, thus upsetting 
tiie cycling procedure. 

The iteration for each plate ir. each aub-cycle, once the CO^ 
array is carpleted. consists of three 3teps. First, a correct weight** 
ang function is determined. This function is designed to penait the 
most rapid change or plate thicknesses in each sub-cycle to eccur on 
platen which will not be changed by the linear extrapolation program. 
These are the plates furthest from the fibers which are found to have 
plates whose primary stress is above the Halting steps. The second 
step is to fill fas second column of that six column array 

described in (CYC. } . This in done in accordance with formula ( i )o 
At the first calling oi this subroutine by (CYC ) the initial thicknesses 
are calculated. This is signified by set ;ing the value of COUNT equal 
to 10. For this case there is no set of thicknesses in column one of 
the array to be averaged, so column two is set directly equal to the 
thicknesses based upon local lateral lead. The last step is the replace- 
ment of column one by the new thicknesses in column two. These are bow 
the new initial thicknesses io be used in the next eub-cycle. Control 
is returned to (CYC.)« 

- 63 » 



LIST OF SYMBOLS GATjCUIATED IK (TBKHES.) 



CO(J) 

‘Hi, CYCIJ3) 

2(1, CYCLE + 1) 

WT 

X 



The plate thickness based upon late- ml load. 
xhe 1hi.cla.ess used fox* each auocyclt of (CYC ^ . 

The thickness found free: the weighted arexage of 
T (1, CYCLE) and CO(J), 

The weighting function for each plate. 

The maximum seam stress associated *ifch each plats. . 
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C 



«- XC i'. .. . „ 

"■y Stfi A ' - .. J '• > ; i ' J Ui i'., . L\. 

dilcKnekses, ^rbltraraly u «•.*.• v,o- *.• - .. a» , ;• y . • J. ,: 

stress; the ecsufctLingrs arc- .♦ t* "i-. •». . iftijiu „ y ^ ^ j. t 

<&ich meets aU c trees crlteri in .. Levin? -.x ds.itgi v.t uu:.- ml. 
cycling procedure is rather isr/clvcd, c so unferstnuri lit desire i^cdo;.;. 
in its entirety, an unierstnn&ir.;; of eneh out -cycle soul its place a 
the overall sche?ae of the min cycle is cauda tory . 

13*8 ultimate alas of each uauin cycle in (CTC . } is to satisix 
the limiting primary beading stress on the* hull girder. Cycle jf-l 
begins with the selection of an arbitrary set of initial thicknesses a* 
mentioned above. This .job is done by subroutine (TliKKSS.) based on a 
stress schedule set up by (STHESI. }«. The stiff ear.* scantlings are now 
determined by (STFN23.) based upon those stresses f on, .1 in (STH5S2. ). 
Subroutine (AM1PI.T.) and (AM1STK) determine the area, moment, aid 
inertia of this initial set of scantlings and now Oycio jfil moves to a 
sub-cycle which designs the scantlings to reals' late: • 1 stress due to 
static water head. {8TKSS1. ) and (STFS32. ) set up a stress schedule 
based upon areas, moments, and inert !:. 0 cf previously calculated refber 
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‘ • N- .iC t 



v iji. : 4 ai 




•L ^ " i • ]L 2 . ..tunic. 






refinemeac; is aeco&pii bed. 




This will be called Cycle 1 , &jr.d reicra to the part of fchs cycle which 
increases the thicknesses of plates that lie on the cuter Tiber. Cycle I, 
will refer to the cycle which attempts to satisfy ths prlwry bending 
stress by increasing plates which would buckle under ths current 
stress schedule. A trigger is set that telio (CYC.) which cycle is 
being executed. Whenever the quanity METE?, :1s set equal tc 1, Cycle A 
is being executed, but when METER equals 0, the subrcu.;ir»te is perform- 
ing Cycle B. Cycles A and 3 both make use of a iincaz program to 
increase the area of the midship cross section to reduce primary bending. 



number of rows equal to the number of plates on the cross-section, The 
first column i3 the thicknesses of plating needed to resist lateral load. 
IThe remaining five columns are set equal to the first, in order that the 
linear program may operate upon the array to chants the thicknesses for 
an optimum stress decrement. 
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L = 1 1>:1 -3 1=' L-") L - 6 




bl 



T. 



tl 



Plates belov neutral axis Increased by lisear 

extragsl&tien . 

Plates abeve neutral axis lac reused by linear 
extrapelatien . 

Caluma index 



FIG. 2 
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*' * ’1 51- . O . . y - i I ( o ]£(■> o..i 

Inc r ernes . ~,t u. ancre-as ■* ' o au ,0 j.l- ;n_r <x »ic oo om ni.’iug in 
co.uuBi 2 anl to the wp plating ». ■ iau r pi-ograra 

calculates the unit clfec of these iKcreiit&jbs anil d a; plays tbe .results 
in columns 4 through c. CtL im is -"he sot cf nt*> cop thicknesses 
causing a stress decrease in the top. Column f> ia the new bottom 
thicknesses causing a stress decrease in the hotter- fibers. Column 6 
has new top and bottoa thicknesses for a dec-ease in. stress at both 
the bottom and top. 

Column 4 through 6 are chec-kel against the enisfeiag stresses 

from Column 1. 



(TOPSTR(l) and ;B0TSTR{i)} to find the ourrect column to use 
as the one vhich lowers the stress for the best attempt to satis!'/ the 

limit stress. 

If Cycle A is being ttua, the top (T. ) or bottom. (T. , } 
sections of the array are plates which are specifically designated . 

They are the bottom plating up to the bilge top, the shear, top and 
stringer atrakes, and the main deck plating. In Cycle 3, tke plates 
arc the top plates which would buckle and tt-.. bottom plates which would 
buckle. During Cycle 3, subroutine (3CK. ) I called to find which 
plates, if any, would buckle for use in the . Inear program. Its secoru 
time it is called to handle e. case whan the Limiting stress is sati"*- 
fled, but scaie plates still, buckle. 
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taoe Gi'cie A uisS y-uc. i fin tu.« iKiur-ncion 

ga:»efe fra*? Cycle A *sua<. be itorc-d v.ai Cycle Bia.ru, Upoa completion 
of Cycle B, "tie r«;suit£ of b ita A u -i £ ar*e jsim,ed ro i thv d«:ii£pae; to 
make a decisis on vhiefc set c :’ tliicJmeeucrs should bs u£t.d as the 
Initial 3i!v in Cycle 2. If ttra dcaiga process is to be tenciaated, 
subroutine (El'll), } is called fo :• final print out of tie dcilga. For 
recycling; all necessary qualities are Initials aed and control is 
transferred back to the refining subcycic for further 1; torsi instability 
calculation . 
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JJC u 


^ £v fc- V. V» 11 . 


hj&'llYZ 


'i’o^aj. rt-nu .',:o 


AR (I) 


uin.-a o. plate * e.s.ii eox'-rr.?. Vj j} 


ASTIF 


Tenporary <.*tor* gr of & rot .a. 


BMfl 


ffe^&uun differs act i>.-o* r e;.~i rn . is: l a.id ac.ua. 
plats fcuc&ling stress i.r cycle i\. 


BDIfS 


Maxiwm. difference bcWMn ro^Okic-d ao.t actual 
plate buShling strrooixi cycle L„ 


BOTSTF(l) 


Beading stre3o on bottaa 1‘tfce tv, ai tea. completion of 
tbs refining sun-eye j.e . 


couiJT 


A count leg device to trigger curtain operations. 


CYCLE 


A count of the desired column ia 


CYCLES 


The current ean cycle of iteration. 


DEF1 


Current deficiency in top streai 


SEES 


Current deficiency in bettor etre s 


DET 


A determinant In this linear prog’vm. 


DT 


Soiall increment ja thickusas oi plates at 
designated locations for linear program. 


XNEET1 


Total tromsnt of Inertia at end at cycle Ac 


JMERT2 


Total iro-tent of inertia at end of cycle B. 


DJRT1A 


Current total moment of ;jiertia of cross section. 


£*(D 


Moment of inertia of each column in T(l.J). 


ISTIF 


Temporary storage of stiffener inertias . 
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ISO IF! 


*<■ i- - ,/ - is r». f) C.G c >ii >; n 


IMD1P2 


Deficiency in priai -5 1c.c .i_g :umm a*; end of • ..» 


METER 


A counting device snowing uivstbcr cycle A or B .'.3 
running. 


MODDK 


Section modulus to deck in cycle a. 


M0D2K1 


Section modulus to deck in cycle b.- 


MODDKL 


Section modulus to ker;X in cycle A. 


MODKL1 


Section modulus to ueel in cycle B. 


MOMEKT 


Current total moment of cross auction. 


1-10(1) 


Total moment of each colurcn in 


MOM1 


Total moment- at end of cycle A. 


M0M2 


Total moment- at end of cycle B» 


S1GS1(I) 


Stress on each plate from cycle B. 


SlGlSV(l) 


Stress on each plate from cycle A. 


SlGL3V(l) 


Initial stress on each plate in each min cycle. 


STOP1 


Trigger to act upon designer's decision. 


TL1M(I) 


The column of T{I,J) used in cycle A. 


TOPSTR(l) 


The bending stress on the top fiber upon completion 
of the refining sub-cycles. 


TSAVE(I) 


The initial thicknesses for cycle* A. 


TS(I) 

x n N 


The initial thicknesses saved for cycle B. 

) 


x i? ; 

x 22 


> Coefficients of the restraint n*t rj s, in the linear 

l program. 
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u(t. ) 



Y2{4) 



n(5) 
rs( 5 ) 



- 1 : ? 

O'.. - 

- iv-l 1" cas.' iiQ , 1 es. c l. ,c. fit. 

is col’.- n .. 

Acru&.l iajj*v*sz is C-. A ’to if s' of t *.• toe fiu*re tr. 
eoluan 5, 

Actual iucreuoe i& oUicJax^ss cf f,e better. HW.rs 
fai co2.\Kin t. 



yi( 6 ) 


Actual increase 
colunn 60 


y 2(6) 


Actual insi’oeae 
An cc,iin iii »>. 



? -1 tulikr as of the top fibers is 
in thiehness of the bo'-t^rn fibers 



I'M. 



The neutral axis of the cross section at cue end 
of each c.vclc. 



o 
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TM« subroutine tents the plat® thicknesses obtained frees 
.local lateral leading for ^stability under the appilcsl ceapreiseive 
stresses . 

The subreut lire la called by (CTC.) and consists; ®£ two 
uai& parts: «ee that tests the longitudinally f rased aecticns for 
buckling and the other that teste the transversely frased section* 
for buckling. 

For sections that era longitudinally fmaed, the* critical 
buckling stress, SIGBUL, is determined far each panel by Bryan 9 s 
Formula. The current panel thicknesses are used in the formal* , sad 
when tw® plate scaas lie betveen longitudinals, tbs subroutlri® 
select® the thickc.es© of the thinner plate to be used. 

The actual primary stress, 31GB1L, is ealculatf/d as the 
stress ©a the point of each panel farthest frees the neutral coda. 

When S1GBKL le less than. S1GB0L, the panel thickness remain* m the 
current value TB. When SiGBKL ie greater than. S1G-MJL, the current 
panel thickness is increased t® the asea® value between the current 
panel thickness, TB, and the required thickness, TL{3) } the thickness 
required to vithstsad SIGBKL® When SISBKL has been decreased to 
vithia 1(W p.s.i. of S13BUL, the panel ia considered to satisfy the 
buckling stress. The factor ef 100 is used te prevent excessive cycling 
ia the emputer. 
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For sections that are transversely framed, the primary stress 
acting es& a panel is calculated toy the subroutine (MAXSTR. ), which 
determines the maximum compressive stress, S1GBXT, acting ©a each 
transversely framed panel. The critical buckling stress, SlGBUY, is 
calculated by Mentgomerie 8 » Formula. S1C-3UT is calculated few each 
plate cetEprisiag the panel, assuming the panel consists only ©f that 
plate thickness • Thus each current plate thickness within the trana« 
verse panel is tested far buckling. If S1GBK9? is less than SlGB?7i ! , the 
plate thicknesses remain as the current thickness T(Z,L). If S1GBKT 
is greater than S1GBUT, the plate thickness is increased to the mean 
value between the current tliickses© T(J,L), nsd the required thickness 
tl(j), the thickness required to withstand S1G3XT . Hfesa S1GBXT is 
within 100 p.s.i. of S1GBUL, the panel satisfies buckling. 

(BUK. ) is executed at three points ia the program: at the 
beginning of Cycle B, and at the end of Cycles A sad B. At the 
beginaiag of Cycle B, (BUX. ) determines which plate panels buckle, so 
that the linear program ia (CYC.) kaows which panel thicknesses to 
change. (BUK.) is executed at the end ©f Cycles A and B t© determine 
if any plates buckle after being changed by the linear program. 

Four possible situations exist when (BUK. ) ie executed at 
the end of Cycle A and B. 

Id Limiting stress satisfied and no plate:? buckle. 

II. Limiting stress Eot satisfied and a© plates buckle. 

III. limiting stress not satisfied and plates buckle. 

IV. Limiting stress satisfied and plates buckle. 



2a case® I and II , the thicknesses catering oo? 

returned to (CYC.). Sr case III, (BUK.) increases the thicknesses ©i 
the plates that buckle 3a the *®naer described previously, and the 
thicknesses are returned to (CYC.), 3 j? ease .XV, (ITJX.) trill cycle 

until mil plates aatisfy the buckling stress. Tfce ihic.'nxe* ses «£ tbs 
pends that buckle are Increased in the mnaer previoualy described. 
Tbs new cr»ss»eectl«aal area properties ©i longitadlnaj. ««t tiauous 
Bsesabers are calculated by (AM&PL5?. ,• ®ri (AM1STP.), u*& the -'res* 
distribution is determined by (3XSSBS1.). Hsw stiffener* qcs-js ngs are 
calculated by (STFS3BR.) end then (TERMS , ) is euc rated t© keep the 
plate thicknesses? up»te~d&te with the new at rosy schedule. The new 
stress distribution and panel thickness eu -Are retarded to the 
beginning ©f (SUK. ), where re-»cycllsg occurs. Only aft*#r til plate 
thlckeacs.es meet the instability requirements will cycling cease, arid 
control returns t© (CYC.). Fur all cases, when the pcs*- 1 thicknesses 
are returned t© (CYC.), they are printed out under the usaci≥ 

" BgKPXHQ GOVERNS *' for Cycle A and 

" BUCKLING QOVgSS for Cycle B. 
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3PXF 

BE 

BOT(M) 

BT 

CEG<I) 

cod) 

S1GBKL 

S1GBUL 

S1CBUT 

TB 

TL{J) 

TOP(K) 

TOICHG 

TP 



LIST OF S SimU vUCJi.: . » 



D.U 2 »nsl«c #f * plate panel la lie?? tr'.'jJi i\e lwi 

Itexlrcua dlffe-reecs bt-tv-s-irt the actual *m 3 critical 
buckling stress acting on. a p* nel* 

B1 Mansion of a plate panel perpendicular le the load* 

The number of a plate wfeica is located below tLr 
neutral axis nnd 'truckles . 

The total number of plates which are located s*elaw t&e 
neutml axis a xtii buckle. 

The number of a plate t bat buckles . 

Tbe cans tent "K r ‘ uaed in Bryan *e Fansula «hica He a 
function of tfc-s edgft restraint conditions and tfce boundary 
conditions . 

Applied primary stream acting on s. langltu&ioolly framed 

panel. 

Critical buckling rtrssi? datejaiiaed from 2r/aa 3 e FoTsjjlau 

Critical bu'kllng strese determined from ^ntgansrle k e 
Formula. 

Actual thickness of a longitudinally framed panel. 

Required t&ickaesa of a longitudinally ©r transversely 
framed panel to withstand - <r uprose tve leading. 

The number of a plate which ia located, above the neutral 
axis and buckles. 

Total number of plates that buckle. 

The total ruober of plate# whicL are located abs»v« the 
neutral axis and buck!.*. 
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